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Understanding the conditions which favor crystallisation or vitrification of liquids has been a long-standing 
scientific problem lOJ-d]. Another connected, and not yet well understood question is the relationship between 
the glassy and the various possible crystalline forms a system may adopt fll^]. In this context, B2O3 is a 
puzzling case of study since i) it is one of the best glass-forming systems despite an apparent lack of low- 
pressure polymorphism ii) it vitrifies in a glassy form abnormally different from the only known crystalline 
phase at ambient pressure |@] iii) it never crystallises from the melt unless pressure is applied, an intriguing 
behaviour known as the crystallisation anomaly 0-E3] • Here, by means of ab-initio calculations, we discover 
the existence of novel B2O3 crystalline polymorphs with structural properties similar to the glass and formation 
energies comparable to the known ambient crystal. The resulting configurational degeneracy drives the system 
vitrification at ambient pressure. The degeneracy is lifted under pressure, unveiling the origin of the crystalli- 
sation anomaly. This work reconciles the behaviour of B2O3 with that from other glassy systems and reaffirms 
the role played by polymorphism in a system's ability to vitrify flcHTill . Some of the predicted crystals are 
cage-like materials entirely made of three-fold rings, opening new perspectives for the synthesis of boron-based 
nanoporous materials. 

PACS numbers: 



Polymorphism, the possibility for a substance to form sev- 
eral distinct crystalline phases of identical composition, is ob- 
served for a wide range of materials. This phenomenon has 
tremendous importance not only per se for understanding the 
crystallisation process but also because of practical implica- 
tions such as the d esig n and control of new materials with 
specific properties 11211 . a major issue for the pharmaceuti- 
cal industry ill 311 . Indeed, the ability of molecular units to 
pack in various ways generates crystal phases which gener- 
ally differ in their physical properties and cohesive energies. 
Another important implication of poly morp hism is related to 
the glassy state: as pointed out earlier OlOi ll ill , glass forma- 



tion is often prevalent for those materials which are found in 
a variety of crystalline forms. An obvious example is silica 
(SiC>2), the archetypal glass-former, which at low pressure is 
found as quartz, cristobalite, keatite, tridymite, coesite and 
moganite IU4I1 . The existence of these many poly types illus- 
trate that the structural units, here the SiC>4 tetrahedra, can 
occur in several conformations with little difference in strain 



energy, allowing for the possibility of metastable states 114 1. 
The ease of glass formation is usually understood as the re- 
sult of the system frustration associated to the presence of 
many minima of comparable energy in the crystal energy land- 
scape (CEL). In the context of organic chemistry, it has also 
been observed that systems with many almost equi-energetic 
structures containing a common interchangeable motif, cor- 
relate with a tendency to disorder 1U5I1 . This is the case 
for instance for Chlorothalonil or Aspirin. On the contrary, 
2,3-dichloronitrobenzene shows a monomorphic behavior be- 
cause of the large energy difference (> 7 kcal.mol -1 ) between 
the most stable form and its other poly types lfl5ll . 



The existence of a correlation between a system's ability to 
amorphise and the presence of numerous crystal polymorphs 
with comparable energies is remarkable since it seemingly ap- 
plies to a very broad range of materials, which includes or- 
ganic, inorganic and metallic systems. In addition, several 
studies have pointed out a structural link between the glassy 
and the known crystalline phases |0,|5t]. These connections 
may also have ramifications for the understanding of phenom- 
ena such as polyamorphism (existence of distinct amorphous 
states), liguid-liquid transitions and even possibly for protein 
folding IU6I1 . 

Boron oxide (B2O3) is one of the simplest glass-forming 
oxides and has received considerable attention Ja, 1 17142 ill . It 



is one of the best glass formers, even a better one than sil- 
ica, since it has never been observed to crystallize from a dry 
melt at ambient pressure: even if the melt is seeded with crys- 
tals and maintained for several months at various temperatures 
below the melting point, no crystal growth is observed at any 
imposed cooling rates (\dT/dt\ > 1CT 5 K.s _1 ) 0-0]. Spon- 
taneous crystallisation from the melt is only obtained when 
the pressure is raised above a threshold level, typically in the 
range 0.4-1.0 GPa JH 0. This puzzling behavior has been 
termed the B2O3 crystallisation anomaly |H|,@|. Although 
nucleation theory based models have been put forward to ac- 
count for this behavior [9], its microscopic origin remains un- 
explained. 

For pressure below 2.0 GPa, the obtained crystal, made of 
ribbons of BO3 planar triangle units, is B2O3-I. At higher 
pressure, another polymorph, B2O3-II, in which the boron 
atoms are four-fold coordinated is crystallised HI. Alterna- 
tively, crystalline B2O3-I can also be prepared by the stepwise 
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FIG. 1: Construction of new B2O3 polymorphs: Vertices of sp - 
carbon structures are replaced either by the BO3 triangle or the 
B3O6 (boroxol) ring as illustrated here in the case of a graphene 
layer. The TO and TO-b polymorphs were obtained by stacking the 
BO3- and BsOs-substituted layers respectively. 



dehydration of orthoboric acid (H3BO3H22] or by seeding a 
melt with borophosphate 12311 . 

The glassy state has also remarkable anomalies: 1) with 
only two experimentally known polytypes, well separated 
in energy, the ease of vitrification of B2O3 seems at odds 
with the correlation mentioned before: a poor polymorphism 
should reflect in a poor glass-former; 2) the glass density ( 1 . 84 
g.cm -3 ) is significantly lower than that of the known crys- 
talline polymorphs, by a magnitude of 30 and 40 % (2.55 
g.cm~ 3 for B2O3-I and 3.11 g.cm" 3 for B 2 3 -H); 3) the 
glass medium-range order is markedly different from that of 
the crystalline polymorphs: in the glass, a majority of the 
BO3 molecular units form three-fold rings, B3O6, known as 
boroxol rings (Fig. 1) llrjl [Til \vh . Strangely, these rings, 
which can be viewed as superstructural units, are totally ab- 
sent in both B2O3 crystalline polymorphs. Boroxol rings are 
however found in large amounts in certain borate crystals such 
as K 3 B 3 6 or Cs 2 0-9B 2 3 JH]. 

The seemingly poor polymorphism of B2O3 compared to 
silica and the fact that pressure is required to crystallize 
B2O3-I from the anhydrous melt raise the question whether 
B2O3-I is indeed the equilibrium crystalline phase at ambient 
pressure. Although unknown polymorphs have been proposed 



in the literature 120112411 . none of them can explain consistently 



the mentioned anomalies since the predicted structures appear 
at too high energies in the CEL. This motivated the present 
systematic exploration of the B2O3 CEL. 

This was done by exploiting the isomorphism between the 
trigonal network of B2O3 and other three-fold coordinated 
networks such as sp 2 -carbon structures (see Methods). We 
investigated all possible networks (13) with up to six vertices 
in the unit cell, as obtained from an exhaustive search origi- 
nally applied to carbon polymorphs [25J|. We also considered 
a layered structure based upon the topology of graphite (Fig. 
1). By placing BO3 triangle units at the vertices of the orig- 




FIG. 2: Examples of nanoporous structures: The B03-made T7 and 
its corresponding B306-made T7-6 structures: large channels, of di- 
ameter approximately 8.0 A and 1 1.0 A respectively, are clearly seen. 

inal structures, 13 topologically different B2O3 polymorphs 
were generated which are labeled TO to T13 in the following. 
Among them, T13 is the known B2O3-I crystal. In T8 and 
T10, 50 % of the boron atoms belong to three-fold rings, i.e. 
boroxol rings. Further to expand the search and to investigate 
the role of the boroxol ring as a structural motif, 1 3 additional 
structures, labeled T0-6 to B2O3-I-6, were generated by re- 
placing the BO3 units by the B3O6 ones. In this way, struc- 
tures which are made of 100 % boroxol units were obtained, 
taking advantage of the self-similarity between a BO3 and a 
B3O6 unit (Fig. 1). All structures were then relaxed using 
first-principles calculations within the density functional the- 
ory framework. 

The resulting B2O3 polymorphs have low (1.0-2.0 g.cm ) 
to very low (< 0.7 g.cm -3 ) densities. Most of the new crys- 
tals are microporous due to cage- or channel-like structures. 
The shape of these channels is either circular, rectangular or 
triangular with a minimum length of typically 6 to 17 A (see 
Fig. 2 for an example). 

Fig. 3a shows the cohesive enthalpy versus density for all 
the obtained crystals. Very strikingly, many of them fall in a 
narrow energy range which includes the known B2O3-I poly- 
morph. This situation, with many minima in the CEL, is ex- 
pected to prevent the nucleation of a given crystal by favoring 
a disordered structure instead. Therefore, this explains the 
glass anomaly 1) (ease of vitrification). 

All new predicted structures have densities lower than 
B2O3-I and the majority of them fall in the range 1.1-1.7 
g.cm -3 , i.e. close to the liquid density above and near the 
glass transition temperature (p/i 9 «id(2000 A') ~ 1.5 g.cm -3 



and , 



s (540 K) ~ 1.8 g.cm 3 ). This is consistent with the 
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FIG. 3: Enthalpy as a function of the density for the obtained B2O3 polymorphs: a) at P= 0.0 GPa b) at P = 1.0 GPa, on the same scale as 
a). The enthalpy reference, Ho, is that obtained for B2O3-I and is shown by the horizontal dashed line. The vertical dashed line indicates the 
glass density. The different symbols and colors refer to the proportion of boroxol rings in the polymorphs. 



observed amorphisation in a low density glass as the temper- 
ature is decreased from the liquid state. This thus explains the 
glass anomaly 2). 

Another result of noteworthy relevance is the fact that the 
majority of the low energy polymorphs are those incorporat- 
ing large amounts (50 or 100 %) of boroxol rings. Indeed, 
the BO3 to B3O6 substitution tend to give polymorphs with 
lower enthalpies (by ~ 2.0 kcal/mol on average, up to <~ 5.5 
kcal/mol in the T3 case) with only two exceptions, T8 and 
T10, which already incorporate 50 % of three-fold rings in the 
starting structures. Interestingly, B2O3-I-6, the boroxol-made 
equivalent of B2O3-I, is found among the lowest energy poly- 
types. These findings, thus, confirm that boroxol rings are 
essential to stabilize low density structures fl9ll and explain 
the glass anomaly 3). We note that in the T10 structure, the 
dihedral angle distribution between adjacent ring and non-ring 
units is well in line with the one measured in the glass by two- 
dimensional NMR lErjll : this supports the idea that the glass 
shares structural similarities with the predicted polymorphs in 
the region of energy-density degeneracy. 

The enthalpy penalty with respect to pressure is larger for 
the new polymorphs than it is for B2O3-I as a consequence 
of their low densities. As a result, the new polymorphs be- 
come less stable than B2O3-I at relatively low pressures (Fig. 
3b). The enthalpy separation between B2O3-I and the sec- 
ond lowest polymorph rapidly increases with pressure (~ 0.8 
kcal/mol at 0.5 GPa and ~ 1.3 kcal/mol at 1.0 GPa), gradu- 
ally leading to a situation in which a monomorphic behavior is 
expected 11511 . This is fully in line with the experimental ob- 
servation that a threshold pressure (0.4-1.0 GPa) is required to 
observe B2O3-I crystallisation from the melt 17H9J]. 

Thus, the present results unravel the mystery of the crys- 
tallisation anomaly: in B2O3, the crystallisation is avoided 
at ambient pressure as a result of the existence of several 
competing phases which eventually induces the system amor- 



phisation. The degeneracy of B2O3-I with the competing 
phases can be lifted by the application of a modest pressure 
which leaves B2O3-I sufficiently separated in energy from 
other competitors. 

These results reveal a much richer polymorphism in B2O3 
than expected before. In this sense, the behavior of B2O3 
appears now to be much more similar to that from other well- 
known glassy systems such as silica. The slight differences 
in energy between the polytypes comes from the variability of 
the B-O-B angle linking the structural units, just as in silica 
Si-O-Si angles bridging the tetrahedra. 

An important issue raised by the present work is whether 
the predicted B2O3 crystals can be experimentally synthe- 
sised and if so, how. All the structures investigated fall within 
the energy range of thermodynamically accessible poly- 
morphs (a debatable quantity, of order of 3 kcal.mol -1 ) 11511 . 
Of course, kinetical or experimental difficulties may however 
prevent their observation. We note that there are at least two 
reports in the literature of low density polymorphs BUI, 
both obtained by dehydration of H3BO3. In the first one 12711 . 
made in the early attempts to obtain B2O3-I, the structure is a 
cubic phase of density 1.805 g.cm -3 . However, in this case, 
the fact that others have not been able to repeat the experi- 
ment remains puzzling. Part of the difficulty may well be due 
to the fact that in this density region, one needs to avoid the 
glass transition resulting from the polymorphic degeneracy in 
the energy-density diagram. The second[28], much more re- 



cent, reports an hexagonal phase of density 0.69 g.cm , in- 
correctly assigned to B2O3-I and which could well be one of 
our very low density polymorphs (< 0.7 g.cm -3 , Fig. 3a). In 
this case, the low dehydration rate used[28] was most proba- 
bly a key ingredient. 

In light of our predictions, further experimental studies ded- 
icated to crystal synthesis would be very much indicated. 
Among the most promising directions towards the synthesis 
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of low-density structures are sol-gel techniques 12911 , Start- 
ing from an alkyl (R) substituted metaboric acid, R3B3O6, 
boroxol-boroxol linkages could be formed via hydrolysis fol- 
lowed by either alcohol or water condensation 12111 . 

Our results confirm that the presence of small rings is a very 
desirable feature to stabilize low density frameworks Il9tl . 
a correlation which was first observed in the context of ze- 
olites lf30h - The propensity of boron atoms to form three- 
membered rings is therefore quite appealing as illustrated by 
the recent discovery of crystalline covalent organic frame- 
works (COF), a novel family of high porosity and low density 
materials 113 111 : COF incorporate small rings such as borox- 
ines, i.e. boroxol rings terminated by a radical such as alkyl, 
alkoxy or triaryl. Thus, our work opens interesting new per- 
spectives for the search and design of boron-based nanoporous 
materials. 

Methods Starting from all possible sp 2 -carbon 3D frame- 
work structures with up to six atoms in the primitive unit 
cell ll25ll . we generated the T1-T14 structures by replacing the 
carbon by boron atoms and placing one oxygen in the middle 
of each B-B bond. The structure numbering corresponds to 
the order of appearance in table 1 of Ref. |25l We rejected 
the structure T12, corresponding to the 6(3)5-09 net in Ref. 
25l since the carbon were almost four-coordinated. The struc- 
ture T14 appeared to be identical to T6 and was not consid- 
ered further. In addition, we constructed a layered structure, 
TO, based on the topology of graphite (Fig. 1). The struc- 
tures TQ-b to T13-6 were obtained from T0-T13 by doubling 
each lattice constant and by substituting each BO3 unit for a 
B3O6 unit. Depending on the polymorph, the primitive cell 
contains from 2 to 27 B2O3 units (10 to 135 atoms). T13 
is the known B2O3-I polymorph, and thus, T13-6 is labeled 
B2O3-T6 in Fig. 3. For each structure, both the atomic po- 
sitions and the lattice cell were optimised using density func- 
tional theory calculations. Preliminary calculations, employ- 
ing molecular dynamics simulations at 500 K were first car- 
ried out with the SIESTA code l32ll . norm-conserving pseu- 



dopotentials and localised basis sets as developed in Ref. |33 
Then, the obtained structures were relaxed at T = K with 



the CASTEP plane-wave code 113411 . ultrasoft pseudopoten- 
tials l35ll and the Perdew-Burke-Ernzerhof generalised gra- 
dient approximation functional ll36tl . We used a plane-wave 
basis-set cutoff of 544 e V and a /c-point grid spacing of 2ir x 
0.05 A -1 or finer. Tolerances on the atomic forces and stress 
were set to 0.005 eV/A and 0.01 GPa, respectively. Because 
of the use of the PBE functional, the equilibrium densities are 
expected to be slightly underestimated. Therefore, the densi- 
ties shown in Fig. 3 have been rescaled by a factor 1.109, as 
given by the ratio of the experimental (2.55 g.cm -3 ) and cal- 
culated (2.30 g.cm -3 ) densities of B2O3-I. For some of the 
networks investigated, several enthalpy minima correspond- 
ing to different symmetries were found. Only the lowest min- 
imum of each topologically unrelated structure is shown in 
Fig. 3. For this reason, we do not show the enthalpy of B2O3- 
0, a recently predicted crystal |20)| of different symmetry but 
the same topology as B2O3-I. Other polymorphs previously 



proposed in the literature appear at higher energies i24ll than 
the upper scale of Fig. 3. 
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